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ABSTRACT 

A  complex  wavelet  of  the  form  V'(*.  y)  =  C{z  —  is  used  in  the  continuous  wavelet  transform  to 

obtain  edges  from  a  digital  image  in  two  orthogonal  directions.  In  the  above  equation  ^(x,  y)  is  a  two  dimensional 
complex  wavelet,  C  and  p  are  constants,  j  is  the  square  root  of  —1  and  x  and  y  are  position  variables.  The  square 
root  of  the  sum  of  the  squares  of  the  real  and  imaginary  parts  of  the  wavelet  transform  are  used  to  find  the  edges  in 
an  arbitrary  direction.  Theoretical  models  are  derived  for  the  cases  of  a  step  edge  and  a  ramp  edge.  The  computer 
implementation  of  the  wavelet  edge  detector  involves  the  creation  of  two  masks;  one  for  the  real  part  and  one  for  the 
imaginary  part.  An  expression  is  derived  for  the  size  of  the  edge  masks  using  the  inflection  points  of  a  curve  derived 
from  the  wavelet.  Finally  the  edge  detector  is  applied  to  two  synthetic  aperture  radar  images  and  the  resulting 
images  are  shown. 


1.  INTRODUCTION 


Edge  detection  in  digital  images  is  a  very  important  operation  in  computer  vision  systems.  There  are  a  large 
number  of  edge  detection  techniques  available,  however  most  of  them  do  not  have  the  capability  of  detecting  edges  at 
different  scales.  Many  of  these  edge  detectors  have  a  small  mask  size  and  do  not  permit  points  on  an  edge  far  away 
from  the  point  under  consideration  to  be  included  in  the  calculations.  Such  edge  detectors  would  include  the  Sobel, 
the  Roberts,  and  the  Laplacian.  These  edge  detectors  also  form  Unite  approximations  of  differentiation  operators 
and  lead  to  results  that  are  sometimes  quite  noisy.  The  purpose  of  this  paper  is  to  present  a  solution  to  the  above 
mentioned  problems  by  use  of  the  continuous  wavelet  transform  (CWT)  in  two  dimensions  which  is  given  below: 


W7(o.,6..a„6,)  - 


(1) 


where  W/(a*,6,,ay,6,)  is  the  wavelet  transform  of  the  function  /(x,y),  rl>  the  wavelet  to  be  chosen 

depending  on  the  problem  to  be  solved,  the  horizontal  bar  over  the  wavelet  in  (1)  indicates  taking  the  complex 
conjugate,  and  Uy  are  the  dilation  parameters,  and  bg  and  5y  are  the  traiislation  parameters.  The  function  /(x,y) 
will  be  used  to  represent  an  image  function  such  that  it  will  be  equal  to  a  specific  graytone  value  at  the  location 
given  by  the  coordinates  (x,y).  It  should  be  noted  that  we  have  taken  a  two  dimensional  function  and  transformed 
it  into  a  four  dimensional  function.  At  first  this  appears  to  have  made  the  problem  much  more  difficult.  The  wavelet 
we  use  in  equation  (1)  is  shown  below: 


V-(x,y)  =  C(x-jy)c-'('’+«'’> 


(2) 


where  x  and  y  are  position  variables,  j  is  the  square  root  of  —1,  and  C  and  p  are  constants  to  be  assigned.  The 
one  dimensional  version  of  this  wavelet  is  equivalent  to  a  first  order  Hermite  function.  The  use  of  Hermite  functions 
as  wavelets  to  be  used  for  edge  detection  was  discussed  by  Grossmann^ ,  but  he  did  not  develop  any  specific  models 
or  provide  any  experimental  results.  This  paper  is  divided  into  six  sections  with  Section  2  discussing  the  analytical 
solution  to  the  step  edge  problem.  Section  3  will  discuss  edge  detection  for  a  ramp  edge.  Section  4  will  derive  the 
wavelet  transform  in  terms  of  the  application  of  two  masks,  one  for  the  real  part  of  the  transform  and  one  for  the 
imaginary  part,  also  an  expression  will  be  derived  for  the  size  of  the  masks.  Section  5  will  describe  some  results  of 
applying  this  wavelet  edge  detector  to  two  synthetic  aperture  radar  images,  and  Section  6  will  present  conclusions. 
The  purpose  of  sections  2  and  3  is  to  clearly  show  the  applicability  of  the  chosen  wavelet  for  edge  detection. 
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2.  SOLUTION  FOR  A  STEP  EDGE 


In  this  section  we  will  derive  the  equation  for  the  wavelet  transform  of  an  image  function  which  consists  of  only 
one  simple  step  edge.  This  step  edge  can  be  represented  by  the  following  expression  for  /(r,  y): 


(A  ifx<0 
=  if.>0 


(3) 


where  A  and  B  are  image  graytones  and  we  take  B  >  A.  The  CWT  and  the  wavelet  rl>{x,y)  were  described  in 
equations  (1)  and  (2).  The  dilation  parameters  a.  and  will  always  be  positive  real  numbers.  When  the  wavelet 
given  by  (2)  is  placed  in  (1)  and  the  CWT  is  separated  into  its  real  and  imaginary  parts  the  following  result  is 
obtained: 


r  r  H..,) 

J-oo  J-oo  L  “X  J 

J— OO*/— OO  L  J 


Putting  the  appropriate  values  for  f(x,y)  into  equation  (4)  will  yield: 

J^ext  oo  \  **»  / 

Now  we  let  «  =  VpK®  ~  ^r)/««]  ^  =  ■v/pKv  “  ^»)/“»]-  Then  (ag/y/p)du  =  dz  and  {ay/y/p)dv  =  dy. 

Making  these  substitutions  in  equation  (5)  we  have: 


C  f' 

Wf{ax,bg,ay,by)  =  / 

<\/W*Uy 


-^e-'di. 

y  '/-oo  vP 


i/  ^e-*dul 

[7-00  P  P  J 

+ ^  r  B J  +  r  ■'4 

7-00  P  [7-00  Vp  >/p  J 

Because  dx  =  0,  the  imaginary  part  of  equation  (6)  becomes  0.  Using  the  following  integrals 

a:e~**(ij:  =  —  ^e”* 

/OO 

e"**  dx  =  y/a 

'OO 


the  real  part  of  equation  (6)  evaluates  to: 

W, 


r/  .  ,  Cyjraxay  (B  -  A)  e“'’(^) 

f{ax,bx,ay)  = - P - 


(4) 


(5) 


(6) 


(7) 

(8) 


(9) 


Equation  (9)  is  the  final  result  for  the  step  edge.  It  is  independent  of  6y  and  when  plotted  as  a  function  of  6c  forms 
a  Gaussian  curve  with  its  maximum  value  centered  at  the  location  of  the  step  edge.  As  the  dilation  parameter 
increases,  the  maximum  value  of  the  Gaussian  curve  increases  and  the  width  of  the  curve  also  increases.  In  other 
words,  as  Oc  increases  the  edge  would  appear  brighter  and  thicker.  The  localization  of  the  step  edge  for  this  wavelet, 
based  on  the  maximum  value  of  equation  (9),  is  exact.  In  general  the  wavelet  transform  will  be  complex  when  the 
step  edge  is  oriented  at  an  arbitrary  angle  with  respect  to  the  x  axis  and  will  consist  of  both  a  real  and  an  imaginary 
part: 

lV/(a*, 6c, flj, 6y )  =  Wfr{ag,bc,ay,b^) jWji{ax,bx,<i]f,bg)  (10) 

The  direction  of  a  unit  normal  to  an  edge  which  points  from  the  smaller  graytone  to  the  larger  graytone  can  be 
computed  by  t^dcing  the  inverse  tangent  of  the  ratio  of  the  imaginary  and  real  parts  of  the  wavelet  transform.  If  we 
let  6  be  the  angle  between  the  x  axis  and  the  unit  normal  to  the  edge,  then  d  can  be  found  as  follows; 


0  —  tan  ^  (Oc  1  I  1  )l^/r  (Oc  i  i  i  ^y  )1  "b 


0*  if  IVyr  is  >0 

180*  if  Wjr  is  <  0 


The  determination  of  the  direction  of  an  edge  from  its  normal  can  be  useful  in  eliminating  edges  in  certain  directions 
or  perhaps  in  emphasizing  edges  in  certain  directions. 

3.  SOLUTION  FOR  A  RAMP  EDGE 

In  this  section  we  derive  an  expression  for  a  ramp  edge.  The  ramp  edge  can  be  represented  mathematically  as 

A  if  I  <  0 

qx  +  b  ii0<x<d  (12) 

B  if  x>d 

The  parameter  d  represents  the  thickness  of  the  transition  zone.  At  z  =  0,  f{x,y)  =  6;  therefore  A  =  b.  At  z  =  d, 
f{x,y)  =  B  =  qd+  A]  therefore  q  =  (B  —  A)/d.  The  dilation  parameters  a*  and  Cy  are  positive  real  numbers.  The 
CWT  will  take  the  form  shown  in  equation  (1).  As  in  section  2  we  let  u  =  \/p(z  —  6,)/ac  and  v  =  \/p(y  —  by)/ay 
making  du  =  {y/p/ax)dx  and  dv  =  i\/p/ay)dy.  The  imaginary  part  of  the  CWT  will  be  equal  to  0  as  it  was  in 
section  2.  The  integral  in  u  can  be  separated  into  three  terms  as  given  below; 

^  C  r  f  AacUC-“’  , 

W /  (<>c  1 6c  I  Oy  1 6y )  —  -  /  dv  \  /  du 

ry/p(d-i.)/a.  r  /q  u  \  1  a-UC-“’ 

+  Br 

P  J 

Using  techniques  similar  to  those  in  Section  2  and  making  use  of  the  definition  of  the  error  function  [erf(z)] 


the  final  form  of  the  CWT  for  a  ramp  edge  becomes: 

=  (15) 

where  h  =  ^[(d  —  6c)/ac]  and  ib  =  —bgy/p/ar.  This  equation  is  independent  of  by  and  when  plotted  as  a  function 
of  6c  will  form  a  curve  where  the  peak  is  located  in  the  center  of  the  linear  transition  zone.  Only  one  edge  is  found 


and  its  brightness  and  width  depends  upon  the  value  for  just  as  it  did  for  the  step  edge  problem  in  Section  2. 
When  Og  increases  the  peak  value  of  the  curve  increases  but  so  does  its  width. 


4.  DERIVATION  OF  THE  WAVELET  TRANSFORM  MASKS 


The  computer  implementation  of  the  wavelet  edge  detector  involves  the  creation  of  two  masks;  one  for  the  real 
part  of  the  wavelet  transform  and  one  for  the  imaginary  part.  The  numerical  values  that  are  placed  in  each  of  the 
masks  can  be  computed  from  the  discrete  approximation  of  the  wavelet  transform,  Wf(ag,bg,ay,by): 


where; 


fm-bg' 

-p 

1  a*  . 

fn-6y‘ 

-pl 

1  ay 

e  * 

It  will  be  assumed  that  each  mask  is  a  square  with  L  rows  and  L  columns.  Requiring  that  L  be  an  odd  number,  we 
let  r  =  m  -  +  1  +  (i  -  l)/2  and  let  s  =  n  -  +  1  +  (L  —  l)/2.  Equation  (16)  can  then  be  rewritten  as; 


Wj {ag,bg ,ay,by)  =  2  /(r  +  6,  -  1  -  (L  -  l)/2,  5  +  6,  -  1  -  (L  -  l)/2)  ■ 

r=l  .=1 

^r-l-(L-l)/2  5-1-(L-1)/2^ 

The  two  masks  Mr(r,s)  and  Af,(r,  s)  can  now  be  written  as: 


-  l-(L-l)/2  s-l-(£-l)/2Y 
o*  ’  /. 


y/a^a^  [  Oj,  J 

In  the  above  two  equations  both  r  and  s  vary  from  1  to  L.  It  can  be  seen  that  the  values  placed  in  the  masks  are 
functions  of  the  dilation  parameters  a*  and  Oy  so  that  the  masks  must  be  recomputed  each  time  there  is  a  change  in 
one  or  both  of  the  dilation  parameters.  The  values  placed  in  the  two  masks  must  also  be  recomputed  if  one  or  both 
of  the  constants  p  and  C  are  changed  or  if  the  size  of  the  masks  is  changed.  If  we  let  Hr  be  the  result  of  applying 
the  mask  Mr(r,s)  to  an  image  at  a  particular  point,  say  {x,y),  and  we  let  Hi  be  the  result  of  applying  the  mask 
Af,(r,s)  at  the  same  point,  then  the  final  computation  for  the  magnitude  of  an  edge  at  this  point  will  be  H[f{x,y)] 
and  is  given  below: 

H[fix,y)]  =  [H^+Hf]^/^  (22) 

The  direction  of  the  normal  to  the  edge  at  the  point  (x,y)  and  measured  with  respect  to  the  i  axis  can  be  represented 
hy  6{x,  y)  and  can  be  computed  as: 


e{x,  y)  =  tan 


0®  if  Hr  is  >0 

180®  if  Hr  is  <0 


The  only  other  thing  which  needs  to  be  derived  is  an  expression  for  the  size  L  of  the  masks  themselves.  This  will 
be  done  by  making  use  of  the  location  of  the  inflection  points  of  a  curve  that  can  be  found  from  the  wavelet,  llie 
wavelet  itself  can  be  divided  into  real  and  imaginary  parts: 


V’(*,y)  =  ^r(*,y)-jV.(*.y)  (24) 

Letting  h(x,y)  =  V*  h(i,i/)  can  also  be  divided  into  real  and  imaginary  parts  [hr(*.y)  and  h,(x,]/)]: 

h(®.  y)  =  Ar(*.y) y)  (25) 

In  solving  for  L  we  will  let  a^,  b^,  y,  ay,  and  by,  be  constants  and  consider  the  curve  which  is  formed  in  the  hr(z,y) 
and  X  plane.  This  curve  will  have  inflection  points,  two  of  which  are  close  to  the  x  axis  and  the  distance  between 
them  can  be  used  in  determining  a  value  for  L.  Letting  u  =  (z  —  bt)/ar  and  v  =  (y  —  by)/ay  and  differentiating 
hr{x,y)  with  respect  to  x  will  yield: 

^  =  -[I  -  2pu*Ie-'’(“"+‘'’>  (26) 

C/X  dg 

Differentiating  hr(z,y)  a  second  time  with  respect  to  z  will  give: 

^  =  -^[2p«  +  (1  -  2pu»)pu]c-'’(“’+“’)  (27) 

An  inflection  point  is  found  by  simply  setting  27  equal  to  zero: 

2pu  +  (I  —  2pti*)pu  =  0  (28) 

One  solution  to  (28)  is  u  =  0.  If  u  is  not  zero  then  the  following  solution  is  obtained: 


Since  u  =  (i  -  6,)/a,  we  have  for  z: 


d^h- 

A  similar  result  can  be  obtained  for  y  by  setting  =  0: 

ay* 


y  =  by  ±ay 


(31) 


The  distance  between  the  two  inflection  points  which  are  the  furthest  away  from  each  other  and  located  along  the 
z-axis  will  be  2a*\/3/(2p).  Along  the  y-axis  the  distance  between  these  inflection  points  will  be  2a,\/3/(2p).  The 
mask  size  must  be  large  enough  to  cover  the  larger  of  these  two  distances  in  order  to  include  the  area  where  the 
wavelet  would  still  have  a  significant  finite  value.  Doubling  the  larger  of  these  two  distances  would  insure  that  the 
masks  cover  an'  area  outside  of  which  the  wavelet  has  an  insignificant  value.  The  value  for  L  must  be  odd  since  it 
will  be  necessary  to  calculate  edges  for  an  image  pixel  at  the  center  of  the  masks.  The  absolute  minimum  mask  size 
will  be  3  X  3  although  in  general  much  larger  mask  sizes  will  yield  better  results.  Given  the  above  constraints,  one 
of  the  following  six  expressions  will  allow  a  value  for  L  to  be  determined  and  the  operator  (int)  means  to  round  off 
to  the  nearest  integer: 

1.  L  =  (int) [da.  \/3/(2p)]  if  a,  >  a,  and  (int)[4arv/3/(2p)]  is  odd  and  >  3 

2.  L  =  (int)[4a,v^3/(2p)]  +  1  if  a,  >  Oy  and  (int)[4ac  \/3/(2p)]  is  even  and  >  2 

3.  L  =  (int)[4ay  v/3/(2p)]  if  Oy  >  a,  and  (int)[4ay  \/3/(2p)]  is  odd  and  >  3 

4.  L  =  (int)[4ay^3/(2p)]  +  1  if  Oy  >  a,  and  (int)[4ay  v^3/(2p)]  is  even  and  >  2 


5.  L  =  3  if  a*  >  fly  and  (int)[4ac  \/3/ (2p)]  <  2 

6.  i  =  3  if  Oy  >  a,  and  (int)[4ay  \/3/(2p)]  <  2. 


5.  APPLICATION  OF  THE  WAVELET  EDGE  DETECTOR 
TO  SYNTHETIC  APERTURE  RADAR  IMAGES 

The  purpose  of  this  section  is  to  show  some  results  of  applying  the  wavelet  edge  detector  to  two  synthetic 
a{>erture  radar  images.  Each  image  consists  of  512  by  512  pixels  and  has  256  gray  levels  (0  to  255).  The  following 
list  of  image  processing  algorithms  was  used  on  the  original  images; 

1.  Speckle  Reduction  using  a  Geometric  Filter 

2.  Edge  Preserving  Smoothing 

3.  Edge  Detection  using  the  Wavelet  Transform 

4.  Thresholding 

5.  Thinning 

6.  Elimination  of  Small  Connected  Components 

7.  Connected  Components 

8.  Outermost  Borders 

9.  Eliminate  Pixels  not  on  the  Outermost  Borders 

When  performing  edge  detection  on  radar  imagery  one  seldom  applies  an  edge  operator  directly  to  the  original 
image.  This  is  because  synthetic  aperture  radar  imagery  contains  speckle  noise  which  is  very  annoying  to  both 
human  interpreters  and  to  computer  vision  algorithms.  A  geometric  filter*  was  used  to  reduce  speckle  noise.  The 
purpose  of  the  edge  preserving  smoothing  algorithm  is  to  eliminate  noise  and  at  the  same  time  preserve  edges 
from  degradation.  The  variation  of  the  gray  tone  in  a  neighborhood  around  each  pixel  is  used  to  determine  the 
direction  that  is  most  homogeneous.  Smoothing  is  then  performed  in  this  direction.  The  specific  scheme  used  for 
edge  preserving  smoothing  was  developed  by  Nagao  and  Matsuyama^.  Eldge  detection  is  performed  using  the  CWT 
as  described  above.  The  following  values  were  used  for  the  constants  C,  p,  o,,  and  Oy.- 

P  =  1 

Os  =  2 

Oy  =  2 

C  =  255(2)p3'V[V  xa^ayw] 

where  w  is  the  smallest  gray  level  difference  for  which  an  edge  will  be  required.  For  the  images  discussed  in  this 
section  w  was  set  equal  to  a  value  of  30.  The  above  expression  for  C  comes  from  the  solution  to  the  step  edge 
problem  which  was  given  in  section  2.  Step  edges  which  have  a  gray  level  difference  of  w  or  greater  will  end  up  with 
the  maximum  gray  value  of  255.  Step  edges  which  have  a  gray  level  difference  less  than  w  will  be  set  to  a  gray  level 
less  than  255.  This  allows  edges  with  small  gray  level  differences  to  easily  be  eliminated  by  a  simple  thresholding 
operation  where  the  threshold  is  set  equal  to  254.  The  value  of  L  (mask  size)  which  corresponds  with  the  above 
parameter  values  is  11.  A  thinning  operation  is  then  required  in  order  to  thin  most  edges  down  to  the  thickness 
of  1  pixel  in  width.  Short  edges  and  components  are  then  eliminated.  The  purpose  of  the  connected  components 
routine  is  to  provide  a  unique  label  for  each  component  of  1-pixels  in  the  binary  image  that  has  been  thinned.  Each 
label  is  a  number  that  is  assigned  to  every  pixel  in  a  given  connected  component.  Use  was  made  of  the  definition  of 
8-connectivity  for  1-pixels  and  4-connectivity  for  0-pixels.  These  definitions  were  also  used  for  the  thinning  operation. 
An  algorithm  which  finds  the  outermost  borders  is  then  utilized.  Border  following  determines  and  uniquely  labels 
all  1-pixels  that  exist  between  a  given  connected  component  of  1-pixek  and  a  connected  component  of  0-pixels.  A 
border  point  is  defined  as  a  1-pixel  that  has  at  least  one  0-pixel  in  its  4-neighborhood.  There  are  two  types  of 
borders,  hole  borders  and  outer  borders.  Outerirwst  borders  are  outer  borders  which  are  not  surrounded  by  any 
other  outer  borders.  The  particular  border  following  routine  used  here  was  developed  by  Suzuki  and  Abe^.  The  last 
image  processing  algorithm  utilized  is  that  of  eliminating  all  of  the  pixels  which  are  not  on  an  outermost  border. 
This  is  a  simple  operation  because  when  a  pixel  is  found  to  be  on  an  outermost  border  it  is  given  a  unique  label 
(255).  Therefore  all  pixels  which  do  not  have  this  label  are  set  to  zero. 


Figure  2:  Radar  Image  of  Airfield  after  Edge  Detection 


Figure  1  shows  a  synthetic  aperture  radar  image  of  an  airfield  located  near  Elizabeth  City,  North  Carolina.  The 
image  is  512  by  512  pixels  and  was  digitized  with  8  bits.  The  radar  system  used  to  obtain  this  image  was  the  UPD-4 
system  which  is  an  X-band  radar  with  1111  polarization.  Figure  2  shows  the  resulting  image  after  edge  detection 
using  the  CWT.  The  edges  defining  the  runways  arc  clearly  visible.  However,  many  edges  located  between  runways 
can  be  seen  which  do  not  correspond  with  any  significant  terrain  feature.  These  edges  can  easily  be  eliminated  by 
the  simple  thresholding  operation  mentioned  above.  Figure  3  shows  the  final  result  after  the  elimination  of  pixels 


Figure  3;  Radar  Image  of  Airfield  After  Elimination  of  Pixels  Not  on  the  Outermost 
Borders 


Figure  Original  Radar  Image  of  a  River  Area 


which  are  not  on  an  outermost  border  The  edges  of  the  runway  pattern  are  clearly  defined  and  have  been  thinned. 
Also  there  are  almost  no  edges  between  runways.  Figure  4  shows  another  synthetic  aperture  radar  image  of  a  river 
located  near  Elizabeth  City,  North  Carolina.  This  image  is  also  512  by  512  pixels  and  was  digitized  with  8  bits. 
The  radar  system  used  to  obtain  this  image  wa.s  also  the  UPD-4.  Figure  5  shows  the  final  result  after  all  nine  image 
processing  algorithms  were  applied  to  the  image  shown  in  Figure  4.  The  edges  defining  the  boundaries  of  the  river 


Figure  5:  Radar  Image  of  River  Area  After  Elimination  of  Pixels  Not  on  the  Out¬ 
ermost  Borders 

are  clearly  seen  while  other  edges  have  been  eliminated. 

6.  CONCLUSION 

The  edge  detector  developed  in  this  paper  using  the  stated  complex  wavelet,  is  a  very  powerful  approach  to  the 
problem  of  edge  detection.  The  theory  for  both  a  simple  step  edge  and  a  ramp  edge  show  that  the  edge  detector  has 
only  one  peak  which  occurs  exactly  where  the  edge  should  be.  Also,  good  results  were  shown  for  some  actual  radar 
images. 
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